Generation of macroscopic singlet states in atomic ensembles 
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We study squeezing of the spin uncertainties by quantum non-demolition (QND) measurement 
in non-polarized spin ensembles. Unlike the case of polarized ensembles, the QND measurements 
can be performed with negligible back-action, which allows, in principle, perfect spin squeezing as 
quantified by [G. Toth et al, Phys. Rev. Lett. 99, 250405 (2007)]. The generated spin states 
approach many-body singlet states, and contain a macroscopic number of entangled particles, even 
when individual spin is large. We introduce the Gaussian treatment of unpolarized spin states and 
use it to estimate the achievable spin squeezing for realistic experimental parameters. Our proposal 
might have applications for magnetometry with a high spatial resolution or quantum memories 
storing information in decoherence free subspaces. 

PACS numbers: 03.65.Ud, 03.67.Mn, 32.80.Qk, 42.50.Dv 



I. INTRODUCTION 

Realization of large coherent quantum systems exhibiting nonclassical behavior is at the center of attention in 
many-body quantum experiments with cold atoms [l[ and ions Q- In a system of spin-i particles spin squeezing 
[H, |i| is one of the most successful approaches for creating large scale quantum entanglement 1 . In a spin squeezed 
state, the collective spin of an ensemble of particles is almost completely set into one direction, while the variance of 
an orthogonal spin component is decreased below the standard quantum limit. In the special case of spin-1/2 atoms 
this implies atom-atom entanglement Q 2 . In a Quantum Non-Demolition (QND) scheme, the particles interact with 
a light field, which is subsequently measured projecting the atoms into a squeezed state [a, H3|- Typically, the length 
of the squeezing dynamics is limited by losses, and a short-time approximation can be applied. In this regime, two 
of the spin components orthogonal to the large mean spin behave like the canonical variables x and p. The Gaussian 



formalism can then be used for modeling 11[ , and can include realistic effects of noise and imperfections [12Hl5l . |17| . 

At this point the questions arise: Is it possible to realize a protocol for the creation of entanglement by squeezing 
the spin uncertainties that would also work for higher spins? This is important since most experiments use atoms with 
larger spins. The solution is not easy: Known methods for creating spin-i entanglement by spin-squeezing cannot 
straightforwardly be generalized to higher spins, without restricting dynamics or the detection to a spin-1/2 subspace 
[Lsl [l9j . Moreover, from the point of view of modeling spin systems, one might ask: Is it possible to extend the 
Gaussian formalism to unpolarized spin states? This is important if we want to depart from the usual setups with 
large collective spin. 

In this paper, we show how to create and detect entanglement by QND measurement of collective spin in an 
unpolarized ensemble. We show how to create a many-body singlet state of atoms without requiring that they interact 
with each other and the system settles in a ground state of some antiferromagnetic Hamiltonian. Our proposal works 
even for particles with a large spin. With realistic experimental parameters for 87 Rb [lol ]. the predicted squeezing 
dynamics are robust to decoherence, and produce a many-atom singlet state. Unlike standard spin squeezing, the 
method creates entanglement even in the limit of very strong interaction, which might be used in experimental 
implementations with cavities (2(i l2l| , or in any multi-atomic system in which a von Neumann measurement of the 
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1 It is instructive to see the relation of spin squeezing to polarization squeezing Q . 

2 For j > 1/2, spin squeezing can be also the result of entanglement between the spin-1/2 constituents of the particles, e.g., entanglement 
between the nuclear and electronic spins Q. 
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collective spin is possible. We demonstrate the validity of the Gaussian approximation for unpolarizcd spin states. 
For the lossless case, we confirm our finding with comparison to the exact model. 

The paper is organized as follows. In Section 2, we present the spin squeezing parameter to detect the entanglement 
of many-body singlet states, and also discuss the properties of the singlets we aim to prepare. In Section 3, we describe 
the squeezing process. First, we consider the lossless case and present a model based on a Gaussian approximation. 
Later, we include decoherence in the model. For the lossless case, we compare our results to the results of the exact 
model. In the Appendix, we present the deltails of the calculations for the exact model. 



II. DETECTING THE ENTANGLEMENT OF SINGLET STATES 



In this paper, we will use the generalized spin squeezing parameter 

f2 (AJ a ) 2 + (AJ,) 2 + (AJ z ) 2 

- j > W 

where J/ are the components of the collective angular momentu m, ( AJ;) 2 = (J 2 ) — (Ji) 2 and for a system of N 



spin-j particles we define J := Nj. It has already been shown in [22H26I ] that any state giving £ s < 1 is entangled 
(i.e., not fully separable). For completeness, we present briefly the proof for (TT]). For pure product states of the form 

l* P ) = ®fcL 1 |V'*) ) wc have 

N 

E (^f- E £(Aj, w )?*>>*;, (2) 

L—x,y,z l—x,y,z k—1 

where denotes the spin coordinates of particle (k) for I = x,y,z. Here we used the fact that ^(Aj^)?^ * > j. 
For a mixture of pure product states, i.e., for separable states, remains true since the variance is concave in the 
state. 

The states giving £ s = are called many-body singlet states [27| . In particular, an equal mixture of all pur e 
singlets, expected to arise in permutationally invariant systems, has intriguing entanglement properties [22l. l24l |25|. 
The bipartite entanglement of this state has already been determined for qubits |28[. It is very mixed, yet its 
entanglement is robust to noise (26j . For qubits, the singlet state studied in this paper is an equal mixture of all states 
composed of two-qubit singlets as can be seen in figure [TJ 

For an imperfect realization, iV£ 2 gives an upper bound on the number of particles uncntangled with other particles 
[13, [H|. This can be seen as follows. Let us consider a pure state of the form t^^L^ipk) <B> \iI>m+i,...,n), which have M 
particles uncntangled with the rest. For such a state, we have 

M N 
(AJ0 2 = £(Ai<% fe> + [A( £ tf% M+1 N) . (3) 

fc=l fc=Af+l 

Hence, 

M 

(AJ x f + (AJ,) 2 + (AJ,) 2 > E( A ^% fc ) + ( A 4%*> + ( A ^% fc > > M 3- (4) 

fe=i 

Due to the concavity of the variance, mixing pure states with M or more uncntangled particles, we still have iV£ 2 > M. 
If N^l for some quantum state is smaller than this bound, then the state cannot be obtained by preparing pure states 
having at least M unentangled spins and mixing them. If JV(1— £ 2 ) is a large number then we can say that entanglement 
between macroscopic number of particles is present in the system in this sense. 

Next, we will review the relevant properties of singlet states. Many-body singlet states are at the center of attention 
in condensed matter physics, and also in other areas of quantum physics. They appear as ground states of many 
fundamental spin models. There have been numerous experiments and experimental proposals for studying such 
states [29l - [3"2l | . Sur pris ingly, the permutationally invariant singlet appears even in quantum gravity calculations of 
black hole entropy |28|]. Its realization is difficult as a T = thermal state of a spin system, since in this case the 
system Hamiltonian is + J 2 + J 2 , which is essentially the sum of two-body interactions connecting all spin pairs. 

Singlets play an important role in quantum information processing. They are invariant under transformations of 
the form Up = exp(— i/3 Jn), which describe the effect of an external homogenous magnetic field on the spins and ft is 
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FIG. 1: For qubits, the permutationally invariant singlet is an equal mixture of all possible arrangements of two-particle 
singlets. Three of such arrangements are shown for eight particles. Note that the eight atoms are arranged in the same way on 
the figures, only the pairings are different. 

the direction of the field. £ s can directly be related to the decrease of fidelity of a pure state |^) under an external 
magnetic field, defined as 

AF:=l-mU p \n 2 - (5) 

Straightforward calculations show that 

AF « (AJ H ) 2 /3 2 < /3 2 Jg (6) 

for small (3 3 . Thus, small £ s indicates insensitivity to external fields, which makes singlet states applicable to encoding 
quantum information in decoherence free subspaces [3~i| , for sending information independent of a reference direction 
[35j . or possibly, for metrological applications in which insensitivity to external homogenous magnetic fields is needed, 
e.g., measuring the spatial gradient or fluctuations of the fields. Note that measuring field gradients is possible with 
two almost c omp letely polarized atomic ensembles in an entangled state [HI, H?} ■ A bipartite singlet can also be used 
for this aim [381 ] - However, using the singlet of a single atomic ensemble can lead to magnetometry with a better 
spatial resolution. 

III. DECREASING THE COLLECTIVE SPIN UNCERTAINTIES 
A. Outline of the squeezing procedure 

In this paper, we consider a procedure to produce an atomic state with £ s < 1 from a non-entangled state with 
£ s ~ 1. For the initial state and throughout the dynamics we have (J) = 0. Suitable initial states include the thermal 
state (completely mixed state), for which £ s = y/j + 1, and pure product states of the form |^ p ) = <S>k\^k), with 
lw/V>*l = J f° r w hich £s = 1- Spin coherent states also have £ s = 1, however, for them (J) y^ 0. After the initial state 
was created, £, is reduced by sequentially measuring and if necessary correcting by feedback the spin components. 

Next, we will examine the dynamics of the expectation values and variances of collective spin operators and we 
will determine how they depend on the particle number N. In particular, A ~ N p will denote that the quantity A is 
proportional to N p for large N. 

The first squeezing step is a QND measurement of the atomic spin component J x . The realization is similar to QND 
measurements in polarized atomic ensembles jjjl IToj|. The atoms interact with light, which is subsequently measured. 
As discussed later, this procedure can reduce the value of (AJ X ) 2 considerably but with only a minor increase in 
(AJ y ) 2 and (AJ Z ) 2 . However, the procedure results in nonzero (J x ) . Fortunately, the expectation value of the spin 
component J x remains small and we obtain | (J x ) | < yN. Feedback, using the data coming from the measurement 
of the light can then be applied to restore (J x ) = 0. This feedback step is similar to the coherent rotation in polarized 
spin-squeezing [lilHtj, but does not require a large average polarization. However, it is now not simply a rotation, 
but an incoherent process. The feedback introduces negligible noise: It can be shown that for making (J x ) zero, 



3 For mixed states, 4J^^ = 4J^ ( (AJ;) 2 bounds from above the Quantum Fisher Information Fq[p, Jfj], which gives the maximum phase 
estimation sensitivity in an interferometer using the quantum state p and unitary dynamics exp(— iJ^t). Here we used the fact that 
F Q [p,Jn]<^Jn? [H. 
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~ VN atoms are affected, introducing extra noise (AJ X ) 2 ~ \/N, which negligible relative to the initial noise in the 
large N limit. 4 . Alternately, instead of feedback, post-selection could be used to identify cases with low (J x ). 

If the initial state was a pure product state of the type mentioned above, measurement of J x is sufficient to produce 
a squeezed state with £ s < 1. For the thermal state, further squeezing steps are needed for J y and J z . Since (J) = 0, 

the generalized uncertainty relations (AJfc) 2 (A J;) 2 > ^- (J m ) do not enforce a measurement back-action. As such, 
the QND measurements can produce a state with (AJ X ) 2 , (AJ y ) 2 , (AJ Z ) 2 all significantly reduced so that £ 2 <C 1. 



B. Lossless case 



We now describe the details of our calculations. We employ methods developed for treating the Gaussian dynamics 
of continuous variable systems [9l-[lll Hij ]. which we adapt to the case of unpolarized ensembles of spin-j particles. 
First, the QND measurement is modeled in the absence of decoherence, which allows comparison with an exact 
calculation. Then, we incorporate decoherence due to light scattering [l3| . 

For modeling the QND pulse, the atoms are described by the J; operators, while the light pulse is characterized by 
the Stokes operators Si [l3|,[l3]- We choose the initial state to be the completely mixed atomic state, go ■= jzfpLW ^ 
and a fully-polarized optical state with (S) = (So, 0, 0). The full system is described by the operators 

D_U Jy. Jz_ _S^_ _Sy_ S* \ (7 \ 

n \ V7' x/J' V7' v^' v^' VsEj {) 

with a covariance matrix 

Trnn ■ — 2 (^rriRn RnRm) (Rm)(Rn)' ($) 

As shown by simple calculations, for large N the initial state is Gaussian for the Rk operators. That is, symmetric 
moments with order higher than second can be obtained from lower order ones according to the theory of Gaussian 



distributions, knowing that cumulants with order three and higher are zero 41(. In other words, concerning the 
moments of Rk, the state is completely characterized by T, (S) and (J). 

The first step of the QND measurement of J x is interaction between the atoms and the light via the Hamiltonian 

H = MIJ X S Z . (9) 

This suggests a characteristic time-scale [H| 

^ : =?^7- ( 10 ) 

The dynamical equations of T mn can be obtained from the Hciscnberg equation of motion for the operators Rk given 
as 

R ( k out) =Rt } -mt\H}, (11) 
with H = 1. For example, the dynamics of R§ is obtained as 

p(out) _ p(in) k p( in ) p( in ) /l n\ 

H 5 - H 5 + H i ■ y lz > 

where the coupling constant is defined as k := ^. Hence, for the dynamics of the variance of i?5 we obtain 

((i? (out) )2) = ((i? (i„) )2) + f i ((i? (in) )2(i? (in) )2) + ^ {R W{ R W,4*> }+ ), (13) 

where {A, B} + is the anticommutator of A and B. Knowing that due to symmetries of the setup for all times (Rk) = 
for k = 1, 2, 3, 5, 6, we obtain 

((AjR (out ))2) = ((Ai? 0") ) 2 ) + «i((Ai?f )) 2 (Ai?[ in) ) 2 ) + ^(R^) 2 ((ARf n) ) 2 ) 



|=(Ai?f n) {A4 in) , A^ in) }+) + ^(AR ( l n) AR^ n) )(R^ ) ) 



(14) 



4 Scattering sets the two-body correlations {j\ j\ ) to zero, where m is one of the spins that went through scattering, while n in one 
of the non-affected spins. From (Jf) > 0, a lower bound on the average two-point correlations is ~ — i. Since ~ v^V atoms are affected 
by scattering, the number of two-point correlations set to zero by scattering is ~ NyN. Thus, the increase of (AJ() 2 is at most ~ \/N . 
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where we used the notation ARk := Rk — (Rk) ■ 

Let us now consider dynamics for t < r. Knowing that (i? 4 ) = -/So for t = and (ARkARiAR m ■ ••) < 1, we can 
examine how the different terms depend on So- We find that on the right hand side of (fT4"| the second term is of order 
■g^, the fourth term is of the order while the rest of the term are of order 1. Thus, assuming a large number of 
photons, that is, a large So, several terms can be neglected and we obtain 

((AjR (out) )2) = ((AjR (in) )2) + g( jR ( i "))2 ((AjR (in))2 ) + 2^ { AR (™) AJ? (in) ) (fl (i„ ) ) _ ^ 

Long, but straightforward calculations show that several terms can also be neglected in the dynamical equation for 
variances of the other Rk variables, assuming large J and Sq. Hence, for the evolution of the covariance matrix 

L P = MT M T (16) 

is obtained, where M is identity matrix, apart from Ms t \ = 4 ^ k. The evolution of the expectation values is 
described as 



Similar analysis shows that the dynamical equations for higher order moments of Rk can also be simplified if t < r. 
When computing the dynamics of these moments, instead of the Heisenberg equation of motion, i?^ out ' = MkiR^ 
can be used. Under this dynamics, which is a linear mapping between operators, the state remains Gaussian. Note, 
however, that this approximation breaks down for much larger times t ~ t^/J. (See also Appendix.) 

The second step of the QND process is measuring S y of the fight. The theory of Gaussian states can be used to 
describe the measurement, which is modeled with a projection [42[ 5 

r M = l p - Tp(p y TpP y ) MP rl. (is) 

Here MP denotes the Moore-Penrose pscudoinverse and P y = diag(0, 0, 0, 0, 1, 0). As described above, the QND 
measurement minimally disturbs the unmeasured components. In the Gaussian approximation the variance of the 
other spin components remain unchanged. 

Measuring the other Jk components is analogous to the J x case. Note that magnetic fields could be used to rotate 
the collective spin to facilitate the measurement of the different spin components. Finally, as we have already noted, 
after each squeezing step, feedback has to be used to restore the zero expectation value of the angular momentum 
coordinates. We have also mentioned that post-selection can be used in the place of feedback. In this case a feedback 
scheme does not have to be realized, however, part of the experiments must be discarded. Thus, when a QND mea- 
surement determines (J;), only cases with | (J,) | < B are retained. We define I(f(x),L) := J_ L f(x) exp(-x 2 /2A 2 )dx, 
where A is the width of the distribution of values obtained when measuring (J;). The effect of post-selection is given 
as var((Jz)) a fter = /zvar(( J/)) before where := [I(x 2 , B)/I(l, B)]/[I(x 2 , oo)/7(l, oo)] and the fraction of post-selected 
experiments is q := 7(1, £?)/I(l, oo). Hence, for a moderate post-selection of B/A = 0.678 we get q = 0.5 for all the 
three squeezing steps, and the variances decrease to 14% of their original values. For the case of B/A = 1.150, we 
get q = 0.75 and the decrease is 37%. 

Let us now make the calculations for realistic parameters. We consider TV = 10 6 87 Rb atoms with spin j = 1, and 
for the light field So = 5 x 10 7 . Sequential squeezing of the x, y and z spin components is shown in figure [2] The 
horizontal axis indicates the total interaction time, with successive intervals of up to 2r for measurement of J x , J yi 
and J z , respectively. Results are shown for squeezing from a thermal state and also from 

Mo--=\+l)^ ®\~J)^- (19) 

We obtain ^ = 0.32 and £| = 0.20 for the completely mixed initial atomic state and for (fTTj|) . respectively. Remarkably, 
the QND interaction can be solved exactly for the initial state ([TTJ]) . as shown in the Appendix. The results are 
presented in figure [2j In that calculation we find that, for large N, time t ~ r x J 25 gives (AJ X ) 2 ~ \fj and the 
two halves of the atoms remain almost fully polarized into the +z and — z directions, respectively. For much longer 
times, the two halves are not fully polarized any more. In particular, for t ~ t-2 := r x \fj we obtain (AJ X ) 2 ~ 1, and 
£f = 4. Thus, we have squeezing even in the long-time (von Neumann) limit. 



5 Equation H18H is analogous to the formula for the behavior of the correlation matrix during a von Neumann measurement in the case 
of Gaussian continuous variable systems described in . However, l|18|l is based on the theory of SU(2) Wigner functions [43[ rather 
than on the theory of Wigner functions for multimode systems. For polarized ensembles, H18H has already been used in [l7ll . 
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FIG. 2: Dynamics of the spin squeezing parameter ^ as a function of the time for 8r Rb atoms with spin 1. Sequential QND 
measurements are made of the x, y and z spin components, with interaction time t up to 2r. The initial state is the completely 
mixed state (upper traces) or the state of (|19[) (lower traces). (Solid) Gaussian approximation without losses. (Dotted) Gaussian 
approximation including losses with (from top to bottom) a = 50, 75, 100. (Dots) Exact model. 



C. Model including losses 



We now incorporate decoherence, following ideas from [l3l. [13. fl7j . adapted to our use of a correlation matrix of all 
the three spin components. In particular, a parameter rj describes the probability that an atom suffers spontaneous 
excitation due to the off-resonant probe, and thus describes the fraction of atoms that decohere during the QND 
process. For simplicity, we assume that the atoms end up in the completely mixed state, which is usual for handling 
the effects of noise. We obtain 



T'p = (1 - r]D)MT M T (t - rjD) + ?y(2 - t))DYv 



(20) 



where D = diag(l, 1, 1, 0, 0, 0) and r no i so = diag(l, 1, 1, 0, 0,0)x j 6 . Here rij is the variance of j x for a spin-j particle 
in a totally mixed state. The decoherence is connected to k through 



K 2 

n = Q-, 

a 



(21) 



where a is the resonant optical depth of the s amp le and Q = 1 for the spin-^ case. For the spin-1 case, Q = | if 
the near-resonant intermediate state has j = [451 ] . Using these techniques, we calculate the degree of squeezing as a 
function of the time for different, experimentally feasible, values of a Il5f. 

The results, shown in figure^ indicate that considerable squeezing is indeed possible under realistic conditions. For 
the completely mixed atomic initial state, we obtained ^ = 0.74,0.61, and 0.54 for a = 50,75 and 100, respectively. 
For comparison, polarized spin states have been squeezed in variance by ~ 50% under similar conditions. Thus, our 
proposal is realistic with present-day technology. The values of a chosen reflect the state-of-the-art for single-pass 
optical probing [l6|, [44| . The use of an optical cavity could boost the effective a by orders of magnitude [2l| . 



IV. CONCLUSIONS 

In summary, spin squeezing of unpolarizcd atomic ensembles by QND measurement shows several intriguing dif- 
ferences relative to polarized samples. The absence of a significant measurement back-action allows simultaneous 



Compare J20} with (13) in fjj 
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squeezing of all spin components and approaches a macroscopic singlet state. Spin squeezing of this type implies 
entanglement of a macroscopic number of particles, for arbitrary spin. To treat this problem, we have extended the 
Gaussian formalism to include the dynamics of all the three spin components and their correlations for an unpolarizcd 
ensemble 7 . The advantage of our approach is that it is possible to determine the area of validity for our model, and 
we can also incorporate von Neumann measurements. In the lossless case, the results agree with an exact calculation. 
Realistic calculations including decoherence indicate that production of these macroscopic singlet states should be 
possible with existing ensemble systems. In the future, it will be interesting to look at the possibility of storing 
quantum information in the decoherence free subspace J47j, formed by different singlet states, obtained for j > h 
when using different initial states for our squeezing procedure. The method of modeling large spin systems can be 
generalized to multiple ensembles and a series of light pulses or even for modeling many-particle quantum systems in 
other areas of physics. 
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Appendix A: Exact model for squeezing (AJ Z ) 2 starting from the state (|19[) 

1. The j = ■x case. 

We briefly describe a method that makes it possible to model exactly the QND process without using the Gaussian 
approximation. First, let us consider the initial state (|19j) for the j = 2 case and measure J x with the QND interaction. 
We divide the atoms into two groups. In the first group, initially N\ atoms are in the | + t-) z state, while in the second 
group N2 atoms arc in the | — 5)2 state. We define the angular momentum operators Jj^j with k = 1, 2 corresponding 
to the two groups. We choose the two halves equal: Ni = N2 = y- Moreover, without the loss of generality we choose 
Nk to be even since in this case Jk,i have integer eigenvalues between — ^ and The initial state can be given in 
the Jk, x basis by 

l*i):=E/iHl m )i,> ( A1 ) 

m 

and 

|*2>:=5>(m)(-ir|m) 2)ai (A2) 

m 

where for large particle numbers f m (x) oc cxp[— ]. Similarly, we define for the state of the light 

:=^.<?(™)MiighM, (A3) 

m 

2 

where g(z) oc exp[— ■§$-]• Hence, we obtain the evolution for the state of atoms and photons as 

|*(t)) = expHJ^t)!*!) ® |* a > ® |*j> 

= EE e "^ 1+ * ) ' n V^' 1 )/ 2 C?'0(- 1 ) , ' a 5(*)|j'i>*lia)x|*)u lt te,». (A4) 
h J2 * 



Modeling lossy dynamics is also possible by writing down the infinite hierarchy of dynamical equations for many-body correlations and 
truncate it at second order [4rj| . 



8 



The projection to the S y = state can be incorporated into the model by introducing w s —: (S z = s\S y = 0). Then, 
the final state of the atoms is 

|*(t)) ex Yl G ^ + 32)fi(h)f2(j2)(-l) j2 \ji) x \j2) x , (A5) 

where G(j) := J] s e"* J j(s)w(s). Here G(j), for large systems, is essentially the Fourier transform of g(s)w(s). 
The value of w(s) matters only for |s| < y/So since for much larger s we have g(s) w 0. For this case, w(s) is to a 
good approximation for successive values of s is alternating between and a constant. Hence, G(j) is very close to a 
Gaussian around j = with a variance ~ J(^) 2 - Thus, t ~ J 25 r gives (AJ^) 2 ~ \/~J, while for £ ~ T2 := rv^7 the 
width of the Gaussian is ~ 1. That is, when computing \^(t)), only states with j\ + j 2 =0 are selected corresponding 
to projecting to the J x = subspace. Thus, for t ~ T2 our setup realizes a von Neumann measurement of J x . For 
large systems the summation in (| A5|) can be replaced by integration. Using these ideas, we obtain 



Jx) « / 



:*(*)lii>xli 



O1+J2) 2 . (A6) 



Due to the absolute value sign, the (— l) 32 term in (|A5[) can be neglected for this calculation, thus the integral of a 
smooth function must be computed numerically for (J 2 ). The dynamics of the other two variances can be computed 
similarly, knowing the matrix elements of M := Jy + J 2 . The change of (M) during this dynamics is negligible. Note 
that [J x , M] = 0, thus a von Neumann measurement of J x does not change (M). 

2. The j > ~ case. 

Finally, the j > | case is analogous to the j = \ case, if we notice that starting from a product state with Ni 
particles in state \j) and N2 particles in state | — j) gives the same dynamics for (Jk.i), as starting from a state with 
2Nij particles in state | + \) and 2N^j particles in state | — \). 
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